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Abstract 

We  present  a unified  approach  to  fast  algorithms  of  various  discrete  trigonometric  trans- 
forms. With  the  help  of  so-called  Euler  formulas  we  describe  an  elegant  and  useful  con- 
nection between  Fourier  matrices  and  trigonometric  matrices.  It  is  known  that  FFTs 
are  closely  related  to  the  factorizations  of  the  unitary  Fourier  matrix  into  a product  of 
uniteiry  sparse  matrices.  Using  these  Euler  formulas  and  FFTs,  we  obtain  fast  algorithms 
of  discrete  trigonometric  transforms.  As  a further  consequence  of  these  Euler  formulas 
and  Gaussian  sums,  we  compute  all  eigenvalues  of  some  trigonometric  matrices. 

1 Introduction. 

The  fast  Fourier  transform  (FFT)  and  related  algorithms  for  orthogonal  trigonometric 
transforms  are  essential  tools  for  practical  computations.  Special  discrete  trigonomet- 
ric transforms  are  the  discrete  Hartley  transforms  (DHT),  discrete  cosine  transforms 
(DCT),  and  the  discrete  sine  transforms  (DST)  of  various  types.  These  transforms  have 
found  important  applications  in  approximation  methods  with  Chebyshev  polynomials, 
quadrature  methods  of  Clenshaw-Curtis  type  (see  [3]),  signal  processing,  and  image 
compression  (see  [4,  6,  9]). 

Euler  formulas  describe  the  algebraic  connection  between  Fourier  matrices  of  a cer- 
tain type  and  corresponding  cosine  and  sine  matrices.  Using  these  formulas,  FFTs  can 
be  transformed  into  fast  and  stable  algorithms  for  the  DCT  and  DST.  Further,  from 
these  Euler  formulas  the  orthogonality  of  various  trigonometric  matrices  follows  imme- 
diately. For  simplicity  we  consider  only  symmetric  trigonometric  matrices,  i.e.  Fourier 
and  Hartley  matrices  of  type  I and  IV  as  well  as  cosine  and  sine  matrices  of  type  I,  IV, 
V and  VHI. 

This  paper  is  organized  as  follows;  first  we  introduce  generalized  Fourier  matrices. 
New  Euler  formulas  for  these  matrices  describe  a close  connection  with  various  ortho- 
gonal Hartley,  cosine  and  sine  matrices.  These  results  simplify  and  extend  former  results 
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of  [9],  pp.  83-96.  Applying  these  Euler  formulas  and  FFTs,  we  obtain  fast  algorithms 
of  discrete  trigonometric  transforms.  As  a further  consequence  of  these  formulas  and 
Gaussian  sums,  we  can  compute  all  eigenvalues  of  orthogonal  symmetric  trigonometric 
matrices. 


2 Euler  formulas  for  Fourier  matrices  of  type  I 

Let  A'  > 2 be  a given  integer.  The  Fourier  matrix  of  type  I is  the  classical  Fourier  matrix 
defined  in  unitary  form 

with  wjv  :=  exp(— 27ri/A).  Note  that  the  Gaussian  sum  (see  [5],  pp.  326-330)  yields  the 
trace  oi  Fh: 
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Glosely  related  with  type  I Fourier  matrices  are  the  cosine  and  sine  matrices  of  types  I 
and  V: 
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Here  we  set  :=  \f2/2  for  f e {0,  A}  and  :=  1 for  j G {1, . . . , A - 1}.  In  this 
notation  a subscript  of  a matrix  denotes  the  order,  while  a superscript  signifies  the  type  of 
the  matrix.  In  the  following.  In  denotes  the  identity  matrix  and  Jjv  the  counteridentity 
matrix,  which  has  the  columns  of  7jv  in  reverse  order.  Blanks  in  a block  matrix  indicate 
blocks  of  zeros.  The  direct  sum  of  matrices  A,  B will  be  denoted  hy  A®  B.  Defining  the 
orthogonal  matrices 
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we  obtain  for  Fourier  matrices  of  type  I the  following  Euler  formulas: 

Theorem  2.1  Depending  on  whether  the  order  of  the  Fourier  matrix  of  type  I is  even 
or  odd,  we  have 
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Proof:  It  is  obvious  that  {I^nY^n  — ^ZN-  Splitting  into  four  blocks 
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and  using  the  classical  Euler  formula  exp(— ia;)  = cos  a;  — isina;,  we  obtain  (2.2)  by 
blockwise  computation  of  {I^nY^n^n-  The  proof  of  (2.3)  is  similar.  □ 

Remark  2.2  An  analogous  result  to  (2.2)  can  be  found  in  [9],  pp.  85-90,  but  with  a 
complex  matrix  instead  of  Compare  also  with  [1].  The  Euler  formula  (2.3)  is  new. 
Note  that  the  results  and  their  proofs  are  simpler  than  in  [9],  pp.  85-90  and  [1]. 
Corollary  2.3  The  matrices  are  orthogonal. 

Proof:  Since  is  unitary  and  is  orthogonal,  C^+i©(— i)Sj^_i  is  unitary  by  (2.2). 
Hence  the  real  matrices  and  Sjj_i  are  orthogonal.  Other  proofs  can  be  found  in 
[4],  pp.  12-16  and  [6]. 

The  proof  for  the  type  V matrices  uses  (2.3)  and  follows  similar  lines.  □ 

Remark  2.4  Results  analogous  to  (2.2)  and  (2.3)  are  true  for  the  Hartley  matrix  of 
type  I (see  [9] , pp.  77-80  and  [8] , pp.  224-227) 
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with  casa:  :=  cos  a;  + sin  a:.  Then  we  obtain  the  formulas 

(^nY^n^n  — ^N+l®^N-l>  (2-4) 

i^^+lY^ZN+l^^+l  — (2-5) 

The  Euler  formula  (2.2)  can  be  used  for  fast  and  numerically  stable  computations  of 
DCTs  and  DSTs  of  type  I:  Let  x € and  y e with  N — 2*  {t  >2)  and  set 


z ^ j e Since  is  real,  we  can  apply  Edson’s  algorithm  for  the  FFT  of 

real  data  (see  [8],  pp.  215-223  and  [7]).  The  output  of  the  conjugate  even  result  is  in  the 
form  U2n^n(^n^)  where  U2N  ■-  Un+i  ® (-i)Fiv-i)  {I2nY-  Therefore  by 
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we  have  calculated  and  Sj^^iy  simultaneously  using  5Nt  flops. 

If  we  have  to  use  an  FFT  with  complex  data,  we  combine  real  data  vectors  x,x'  G 

(X  I \ 

y + iy'  )'  Then  we  can 

compute  two  DCTs  Clf^iX,Cjlf_^ix'  and  two  DSTs  S)v_iy,  simultaneously  via 

an  FFT  of  length  2N  applied  to  the  complex  input  vector 

In  a similar  way,  the  Euler  formula  (2.3)  can  be  used  for  fast  computations  of  DCTs 
and  DSTs  of  type  V:  For  given  a:,®'  € R^+*  and  y, y'  e R^  the  transformed  vectors 
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, S]^y,  and  S]^y'  can  be  calculated  at  the  same  time  as  components  of 
^N+i^^+i^'  = (^+1  ® where  we  use  an  FFT  of  length  2N  + 1 

with  complex  data  I^+iz'.  If  2N  + 1 = 3*  or  more  generally,  if  2N  + 1 is  a product  of 
small  primes  (see  [8],  pp.  76-101  and  [7])  the  FFT  of  length  2N  + 1 can  be  computed 
very  efficiently. 

3 Euler  formulas  for  Fourier  matrices  of  type  IV 

The  Fourier  matrix  of  type  IV,  defined  by 
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is  related  to  the  Fourier  matrix  of  type  I by  the  formula 
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with  Wn  :=  diag(o;2jv)^o^  i®  therefore  unitary.  If  V is  a power  of  2 or  3,  then 
can  be  factorized  into  a product  of  sparse  unitary  matrices. 

Lemma  3.1  The  trace  of  the  Fourier  matrix  of  type  IV  is  equal  to 
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Proof;  We  begin  with  the  generalized  Gaussian  sum  (see  [5],  p.  330) 
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which  we  split  into  two  sums  containing  even  and  odd  j respectively.  Then 
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and  the  results  follows  by  (2.1).  □ 

Now  we  introduce  cosine  and  sine  matrices  of  type  IV  and  VIII  which  are  closely 
related  with  the  Fourier  matrix  of  type  IV; 
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As  above  we  define  orthogonal  matrices 
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Theorem  3.2  For  the  Fourier  matrix  of  type  IV  and  even  resp.  odd  order,  we  obtain 
the  following  Euler  formulas: 


(3.3) 

= qTe(-i)sr+i-  (3.4) 

Proof:  Similar  to  that  of  Theorem  2.1.  □ 

Corollary  3.3  The  matrices  ,S^ ,C^  and  are  orthogonal. 

Remark  3.4  An  analogous  result  to  (3.3)  can  be  found  in  [9],  pp.  94-96.  Compare  also 
with  [1].  Formula  (3.4)  is  new.  A different  proof  of  the  orthogonality  of  and  can 
be  found  in  [6]. 

Remark  3.5  Similar  formulas  as  in  Theorem  3.2  are  true  for  the  Hartley  matrix  of 
type  IV  (see  [1,  2]) 
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The  Euler  formulas  can  be  used  for  a fast  and  numerically  stable  computation  of 
DOT  and  DST  of  types  IV  and  VIII: 

Using  (3.3)  and  (3.1),  for  arbitrary  x,x',y,y'  € the  DCTs  C^x,  C^x'  and 
DSTs  y and  y'  can  be  calculated  via  one  FFT  of  length  27V  with  complex  data 

I^z'  and  z'  ^ ^ j.  If  TV  = 2*,  this  procedure  requires  about  lOiVt  operations. 

Likewise  by  (3.4),  for  x,  x'  e ,y,y'  € R^+i  the  DCTs  of  type  VIII,  C^x,  C^x' 
and  the  DSTs  S]^iy,  S^iy'  can  be  calculated  via  one  FFT  of  length  27V  + 1 with 
complex  data  I^+iz'. 

Remark  3.6  The  sine,  cosine.  Hartley,  and  Fourier  matrices  considered  above  enjoy 


the  interesting  intertwining 

relations  (see  [2]): 
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with  the  diagonal  matrix  Sw+i  :=  diag((-l)'')f^o  reflection  matrix  J)v  “ 

1©  Jiv-i.  Therefore  applying  (3.7)  in  the  above  algorithm,  it  is  also  possible  to  compute 
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four  DCTs  (or  four  DSTs)  of  type  IV  and  order  N via  one  FFT  of  length  2N  with 
complex  data. 

4 Eigenvalues  of  trigonometric  matrices 

Finally  we  determine  the  eigenvalues  of  trigonometric  matrices  introduced  above.  Since 
the  cosine  and  sine  matrices  of  type  I,  IV,  V and  VIII,  and  the  Hartley  matrices  of  type 
I and  IV  are  real,  symmetric  and  orthogonal,  only  1 and  -1  are  possible  eigenvalues. 
For  a;  e M we  denote  by  [a:J  fesp.  [x]  the  integer  fc  6 Z with  k < x < k + 1 resp. 
k — 1 < X < k. 

Theorem  4.1  The  sine  and  cosine  matrices  , C^,  and  of 

order  N >2  possess  the  eigenvalues  1 and  —1  rvith  multiplicities 

m{l)=\N/2],  m(-l)=  [V/2J. 


Proof:  Since  is  symmetric  and  orthogonal,  only  1 and  —1  can  be  eigenvalues.  Their 
multiplicities  fulfil 

m{l)  + m{—l)  = N. 

On  the  other  hand,  since  clf  and  Sl;_2  are  real,  it  follows  from  (2.2)  and  the  trace 
formula  (2.1)  that 


m(l)  - m(-l)  = trC^  = Re(tri^^_2)  = Re 


l + i2^-2 
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1 for  odd  N, 
0 for  even  N. 


From  these  two  linear  equations  we  obtain  m(l)  = \N/2]  and  rn{—l)  = [iV/2j.  In  the 
other  cases,  the  proof  is  similar.  □ 


From  Theorem  4.1  and  the  Euler  formulas  (2.2)-(2.3)  and  (3.3)-(3.4)  it  follows  im- 
mediately: 


Corollary  4.2  The  Fourier  matrices  of  type  I and  IV  have  only  eigenvalues  1,  —1,  i,  — i 
with  multiplicities: 
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From  Theorem  4.1  and  formulas  (2.4)-(2.5)  and  (3.5)-(3.6)  it  follows: 

Corollary  4.3  The  Hartley  matrices  of  type  1 and  TV  have  only  eigenvalues  1 and  — 1 
with  the  following  multiplicities: 
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